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Abstract. We report the energy dependence of normal branch oscillations (NBOs) in Scorpius X-1, a low-mass X-ray binary 
(LMXB) Z-source. Three characteristic quantities (centroid frequency, quality factor and fractional root-mean-squared (rms) 
amplitude) of a QPO signal as functions of photon energy are investigated. We found that the NBO centroid frequency decreases 
with increasing photon energy when it is below 6-8 keV and turns positively correlated with photon energy at the higher-energy 
side. The rms amplitude increases significantly with the photon energy below 13 keV, and then flats out in the energy band of 
13-20 keV. There is no clear dependence on photon energy for the quality factor. Based on these results, we suggest that the 
NBO originates mainly from the transition layer. 
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1. Introduction 

The bright persistent neutron star low-mass X-ray binary 
(LMXB) Scorpius X-1 traces a "Z" track in the X-ray color- 
color diagram (CCD), which consists of three branches — the 
horizontal branch (HB), the normal branch (NB), and the flar- 
ing branch (FB, Hasinger & van der Klis 1989). The typical 
Fourier power spectra in Scorpius X-1 are composed of distin- 
guishably aperiodic variabilities in most of these branches, in- 
cluding noise components (broad structure) and quasi-periodic 
oscillations (QPOs, narrow feature). There are three distinct 
types of QPOs observed by Rossi X-Ray Timing Explorer 
(RXTE), i.e., the normal branch oscillation (NBO) at ~ 6 Hz, 
the horizontal branch oscillation (HBO) at ~ 45 Hz with a har- 
monic at about 90 Hz, and the kiloherz (kHz) QPOs (van der 
Klis et al. 1996, 1997). 

The NBOs in Scorpius X-1, which occur on the mid- and 
lower NB, show some correlations with the variability of the 
~ 45 Hz HBO and the twin kHz QPOs, (van der Klis et al. 
1996; Yu 2007), and therefore imply some kind of coupling be- 
tween the three types of QPOs (van der Klis et al. 1996; Dieters 
& van der Klis 2000). The NBO frequency remains approxi- 
mately constant at about 6 Hz, in general. However, the timing 
properties along the Z track shows that the NBO frequency ex- 
tends to at most ~ 21 Hz and its position moves smoothly into 
the lower part of the FB (Priedhorsky et al. 1986; Dieters & 
van der Klis 2000). This smooth transition indicates that the 
NBO and FBO are physically related to each other (Casella 
et al. 2006). Properties of both the twin kHz QPO and 45 Hz 



HBO depend on the NBO flux. The upper kHz QPO frequency 
is anticorrelated with the NBO flux, and the lower kHz QPO 
becomes stronger when the NBO flux is low (Yu et al. 2001). 
Significant HBOs are detected during the NBO phase of high 
flux, while the HBO disappears at the low flux of NBO phase 
(Yu 2007). The coupling between the properties of kHz QPO 
and HBOs and the phase of the NBO makes the NBO a unique 
phenomenon. 

The photons from different regions of the accretion disk 
carry different energies and present distinct physical properties. 
So the QPOs detected in different energy band may have vary- 
ing characteristics. Inspired by the idea that the photon energy 
dependence of QPOs may provide some additional information 
for the nature of QPOs, and limited by the available data which 
provide enough spectral information and at the same time in- 
clude QPO signals, we perform a systematically investigation 
for the energy dependence of NBOs in this work. We describe 
our data reduction and analysis in section 2. Results are present 
in section 3. In section 4, we discuss physical implications of 
our result. Section 5 contains a short summary. 

2. Data reduction and analysis 

To investigate the energy dependence of QPOs in Scorpius X-1, 
high time resolution data with enough spectral information are 
needed. QPO signals should also be present in these data. We 
searched all the RXTE observations of Scorpius X-1 for suit- 
able data. Among all the RXTE Sco X-1 data, two binned data 
modes, i.e., B.4ms_16A_0.249_H and B.4ms_16B_0.249_H, 
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which have a 4-ms time resolution and 16 energy bands com- 
posed of all the PCA channels (from original channel to 249) 
binned with two different schemes (A and B), provide adequate 
timing and spectral information for our purpose. We examined 
all the data recorded in those two modes and selected those hav- 
ing QPO signals present (with quality factor Q larger than 2) in 
the energy band of channel 0-63, which roughly corresponds to 
2 to 20 keV. Those selected data are listed in Table [1] We fur- 
ther defined 6 energy bands, as shown in Table|2] for the energy 
dependence study. An HBO signal was found in the highest en- 
ergy band defined in Table in the data set 30035-01-07-00, 
but cannot be detected in the original 0-63 channel. In addi- 
tion, another ~ 45 Hz HBO signal was detected in the original 
0-63 channel in 93082-01 -02-05, but cannot be found in any in- 
dividual energy band defined in Table[2] We also found another 
data set, 30035-01-01-00, having the ~ 45 Hz HBO signal in 
the original 0-63 channels but without NBO. In futher exam- 
ination on this data set, its HBO signal was only detected in 
the original 0-16 channels but not in higher energy bands. We 
therefore did not include this data set in this work. Limited by 
the data available, we study only NBOs in this work. Some of 
those data sets include several segments. We only exploited the 
segments which show the NBO signal (see TableQ]for details). 

The data were reduced with the HEASOFT package ver- 
sion 6.8. Only those data with the elevation angle > 10° and 
the time since South Atlantic Anomaly (SAA) < or > 30 min 
were selected for further analysis. We extracted light curves 
with a time resolution of 4 ms (2~ 8 s) in the six energy bands 
defined in Table|2] Using the Leahy normalization (Leahy et al. 
1983), we computed the PDS without subtracting white noise. 
The average power expected from a Poisson distribution is 2. 
The original 4 ms (2~ 8 s) time resolution was used, and each 
continuous set was divided into segments of 8192 bins. A fast 
Fourier transform (FFT) algorithm for each segment was com- 
puted yielding the frequency PDS in the range of 3 x 10~ 2 - 128 
Hz. We calculated the average PDS for each continuous data in- 
terval. A logarithmic rebinning of -1.02 was used for the rebin 
option. Figure[T]shows an example of the ~ 6 Hz NBO signal. 

We employed three components to fit the resulting PDS: a 
Lorenzian for the very low frequency noise (VLFN), a power 
law for the white noise, and another Lorentzian for the NBO. 
All the fits were performed using the x 2 minimization tech- 
nique. The errors of the analytic model parameters were ob- 
tained at the^f 2 variation of 2.70, which corresponds to a 90% 
confidence for a single parameter. 

A QPO signal is characterized by three quantities: the cen- 
troid frequency (representing the position of a QPO), the qual- 
ity factor (characterizing the coherence of a QPO signal), and 
the fractional root-mean-squared (rms) amplitude (a measure 
of the signal strength, which is proportional to the square root 
of the peak power contribution to the PDS). From the fitting 
we obtained three parameters of NBOs, i.e., the oscillation 
amplitude, the centroid frequency, and the full width at half 
maximum (FWHM). The quality factor was then computed 
according to its definition (centroid frequency/FWHM). The 
rms amplitude was calculated by the integral of the normalized 
Lorenzian PDS in the appropriate frequency range and taking 
into account the mean count rate of the source. Details of rms 



Table 2. RXTE/PCA energy bands used in this work. 



The observations in 1996 and 1998 



PCA 


Energy Band 


Centroid Energy 


Channel 


(keV) 


(keV) 


0-13 


1.83-5.07 


2.54 


14-19 


5.07-7.24 


6.16 


20-22 


7.24-8.34 


7.79 


23-28 


8.34-10.52 


9.43 


29-35 


10.52-13.09 


11.81 


36-63 


13.09-23.46 


18.26 


The observations in 2009 


PCA 


Energy Band 


Centroid Energy 


Channel 


(keV) 


(keV) 


0-11 


2.18-5.26 


2.63 


12-15 


5.26-7.02 


6.14 


16-19 


7.02-8.78 


7.90 


20-22 


8.78-10.10 


9.44 


23-30 


10.10-13.60 


11.85 


31-63 


13.60-28.40 


21.00 



93082-01-01-03 



Frequency (Hz) 

Fig. 1. One example of PDS which shows a ~ 6 Hz NBO. The 
PDS is produced in the PCA original channels 0-63. The mean 
count rate is 14348 cts/s/PCU. The exposure time is 938 s. 

determination can be found in Lewin et al. (1988) and van der 
Klis (1989).. The errors of the quality factor and rms amplitude 
were estimated with standard error propagation (Bevington & 
Robinson 2003). 

3. Results 

We investigated the photon energy dependence of three char- 
acteristic quantities (i.e. the centroid frequency, quality factor, 
and rms amplitude) of NBOs. Figure [2] shows the energy de- 
pendence of the centroid frequency. Although error bars are 
large, there seems to be a common trait in almost all the data 
sets that the dependence evolves from a negative correlation 
to a positive one with increasing photon energy. We therefore 
adopt a two-power-law function, v(E) = a E r ' + b E Y2 , to de- 
scribe the relation between the centroid frequency v (in units of 
Hz) and the photon energy E (in units of keV). The fitting re- 
sults are listed in Table [3] In all the data set we have a negative 



J. Wang et al.: Energy dependence of NBO 



3 



Table 1. Data log. Here listed are the RXTE observation identification number (ObsID), their start and stop time, the exposure 
time in seconds of the light curve after screening (on elevation and time_since_SAA) and selection of NBO signal presence, the 
number of working PCU, the original 0-63 channels count rate in units of cts/s/PCU and the segments with NBO signals. 



ObsID 


Start 


End 


Exposure 


PCU 


Rate 


Seg 








(s) 




(cts/s/rL,U) 




10056-01-04-00 


27/05/96 20:42:05 


27/05/96 23:02:13 


2633 


5 


18776 


All 


i nn^A m c\a iy) 
1 UU JO-U 1 -U4-U2 


zo/uj/yo u4.zo.zj 


ZS/Uj/yO LU.Zy.Lj 


llOUU 


c 

J 


1 / Do J 


A 11 
All 


10056-01-05-00 


28/05/96 11:11:04 


28/05/96 14:23:13 


5277 


4 


16469 


All 








375 


5 






30035-01-02-000 


31/05/98 00:26:47 


31/05/98 08:26:47 


77 


3 


18182 


4,9 








2936 


5 






30035-01-05-00 


01/06/98 00:26:55 


01/06/98 01:18:14 


656 


3 


18413 


All 


30035-01-07-00 


02/07/98 06:26:32 


02/07/98 10:27:14 


4968 


5 


16486 


1,2 


93082-01-01-01 


16/06/09 23:43:14 


17/06/09 00:56:18 


48 


4 


14718 


All 








926 


5 






93082-01-01-03 


18/06/09 01:40:22 


18/06/09 02:21:18 


938 


5 


14348 


All 


93082-01-01-04 


17/06/09 05:13:12 


17/06/09 06:01:18 


1556 


5 


14486 


All 


93082-01-02-05 


23/06/09 04:07:09 


23/06/09 04:40:18 


1311 


3 


13447 


All 



and a positive power index, which together make the centroid 
frequency decrease with increasing energy first and then when 
beyond about 6-8 keV it turns to be positively correlated with 
the photon energy. This dependence is not obvious in one data 
set (ObsID 30035-01-05-00) among the ten under study. While 
the small xl value indicates that this description is only sug- 
gestive due to the large error bars, the appearance of this trait 
in almost all the data sets asserts its reality. Figure [3] shows 
the energy dependence of the quality factor. We do not see any 
clear relation between the quality factor and the photon energy. 

In Figure |4] we can see that the rms amplitude increases 
monotonically with the photon energy except for the highest 
energy band. It means that the strength of the NBO signal in- 
creases with the photon energy. The rms amplitude increases 
significantly with the photon energy below 13 keV. However, 
the increase of rms amplitude with photon energy seems to 
stop when the photon energy is higher than ~ 13 keV, and there 
might be a plateau. We do not have enough photons at even 
higher energies for our study. It would be interesting to know 
beyond which energy the rms amplitude starts to decrease. To 
describe the increase of the rms amplitude with photon energy, 
we use a power law relation, rms(£) = C E r , to fit the data 
points in the lower five energy bands. Table|4]shows the fitting 
results. Roughly speaking, the rms amplitude increases with 
photon energy according to a power law of power index about 
0.7 before reaching a possible plateau. In the data set 30035- 
01-07-00 the plateau seems to start at a lower energy. 

4. Discussions 

The mechanism of NBOs was suggested to be related to the ra- 
dial oscillation during radiation feedback (Fortner et al. 1989; 
Miller & Park 1995) or the different modes of disk oscillations 
(Alpar et al. 1992; Titarchuk et al. 2001; Bildsten & Cutler 
1995; Bildsten et al. 1996; Horak et al. 2004). The standard 
explanation attributes the ~ 6 Hz NBO to part of the accretion, 
with near Eddington accretion rate, in an approximately spher- 



ically symmetric radial inflow at a radius of ~ 100 km (van 
der Klis 1995; Hasinger et al. 1990). According to the radia- 
tion hydrodynamic model (Fortner et al. 1989), NBO originates 
from a radiation-force/opacity feedback loop within a spheri- 
cal flow region at about 300 km from the neutron star surface. 
The frequency of the oscillations depends on the luminosity. 
Alternatives are the acoustic oscillation in a thick disk due to 
the density and optical-thickness perturbations caused by the 
rotating medium in a subsonic region of the accretion disk (the 
frequency is determined by the rotation rate and vorticity, Alpar 
et al. 1992), g-mode oscillations arising from the thermal buoy- 
ancy in the ocean on rotating neutron stars (Bildsten & Cutler 
1995; Bildsten et al. 1996), the low-frequency modulation in a 
non-linear resonance oscillation of the relativistic disk (Horak 
et al. 2004), and the acoustic oscillations of a spherical shell 
surface within ~ 20 km of the neutron star surface because of 
the change of accretion geometry in the transition zone at high 
accretion rate (Titarchuk et al. 2001). All of these interpreta- 
tions focus on the ~ 6 Hz frequency and its dependence on 
the accretion rate. None of them can explain the energy depen- 
dence of the centroid frequency and rms amplitude. 

In the NBOs of Scorpius X-l, both the centroid frequency 
and rms amplitude display strong energy dependence. The cen- 
troid frequency changes from anticorrelation to a positive cor- 
relation at about 6-8 keV. Such a nonmonotonic energy depen- 
dence of the centroid frequency implies that the emission zone 
of the NBOs experience a radial variation during the accre- 
tion process. The rms amplitude increases monotonically and 
significantly with the photon energy below ~ 13 keV. When 
the photon energy is higher than 13 keV, the rms amplitude 
seems to stop increasing. It may indicate that the NBO sig- 
nals reach the strongest strength at 13-20 keV. Higher energy 
photons are often emitted from the inner region of accretion 
disk. Consequently, the most possible region responsible for 
such physics can be referred to the transition zone between 
the inner boundary of accretion disk and the magnetosphere, in 
which the transition of the properties for accretion flow occurs 
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Table 3. Fitting results of a two-power-law description (v{E) = a E r ' + b E Fl ) for the energy dependence of the centroid frequency 
of NBOs in Scorpius X-l. 



ObsID 




a 




r 2 


b 


5 

Xv 


10056-01-04-00 


-0.17+0.058 


3.40= 


=0.65 


0.11 ±0.028 


3.05±0.39 


0.19 


10056-01-04-02 


-0.22+0.026 


3.65= 


=0.29 


0.15+0.030 


2.79±0.43 


0.13 


10056-01-05-00 


-0.62+0.023 


2.52= 


=0.30 


0.086±0.025 


5.03±0.38 


1.03 


30035-01-02-000 


-0.42±0.038 


2.32= 


=0.43 


0.083+0.022 


4.38±0.29 


0.02 


30035-01-05-00 


-0.43+0.11 


0.81= 


= 1.13 


0.040+0.054 


5.23±0.77 


0.65 


30035-01-07-00 


-3.58+0.27 


4.27= 


=0.30 


0.016±0.014 


6.22±0.20 


0.37 


93082-01-01-01 


-0.24+0.052 


5.87= 


=0.70 


0.35±0.050 


1.38±0.68 


0.36 


93082-01-01-03 


-0.16±0.084 


3.16= 


= 1.00 


0.082+0.029 


3.65±0.46 


0.15 


93082-01-01-04 


-0.41+0.057 


3.52= 


=0.85 


0.10±0.042 


4.19±0.68 


0.62 


93082-01-02-05 


-0.29+0.067 


3.11= 


=0.78 


0.11 ±0.025 


3.56±0.35 


0.29 



Table 4. Fitting results of a power law (rms(£) = C E r ) for the increase of rms amplitude of NBOs in Scorpius X-l. Only the 
five lower energy bands are included in the fitting. 



ObsID 


r 


C 




10056-01-04-00 


0.78±0.18 


1.06±0.40 


0.27 


10056-01-04-02 


0.71±0.09 


1.92±0.34 


0.44 


10056-01-05-00 


0.75±0.06 


1.77±0.22 


0.76 


30035-01-02-000 


0.73±0.12 


1.27±0.30 


0.30 


30035-01-05-00 


0.63±0.34 


1.12+0.80 


0.14 


30035-01-07-00 


0.70±0.07 


1.79±0.25 


3.13 


93082-01-01-01 


0.70±0.19 


1.52+0.61 


0.36 


93082-01-01-03 


0.74±0.20 


1.52±0.66 


0.28 


93082-01-01-04 


0.72±0.18 


1.44±0.56 


0.18 


93082-01-02-05 


0.70±0.16 


1.65±0.55 


0.24 



(Titarchuket al. 1998; Titarchuk & Osherovich 1999; Titarchuk 
et al. 1999). NBO occurs on high accretion rates, i.e. near the 
Eddington rate. The high mass flux makes more and more mat- 
ter deposit in the transition zone, leading to the expansion and 
the radial scale change of this region. However, the pile of more 
deposit with greater viscosity can contribute to the radiation en- 
hancement, as well as the emission of photons with higher en- 
ergies. As a result, the radial scale change of the transition zone 
may be responsible for the observed nonmonotonic energy de- 
pendence of the centroid frequency. In addition, because of the 
limited region of the transition zone, more deposit enhances 
the viscosity and thus suppresses the further monotonicly in- 
creasing strength of NBO signal, which can explain the energy 
dependence of rms amplitude at high energy. We therefore sug- 
gest that the NBOs be one kind of oscillations in the transition 
zone. 



5. Summary 

We investigated the energy dependence of three characteristic 
quantities of NBOs in Scorpius X-l. We have the following 
conclusion: 

(1). The energy dependence of centroid frequency changes 
from a negative correlation to a positive one. The turning points 
are located at 6-8 keV. 



(2) . It seems that the quality factors have no obvious rela- 
tion with the photon energy, which indicates that the coherence 
of NBOs is independent on the photon energy. 

(3) . The rms amplitude increases significantly with the pho- 
ton energy below 13 keV. It seems to reach a plateau beyond 
that. 

(4) . The observed energy dependence of three characteris- 
tic quantities of NBOs in Scorpius X-l implies that the NBOs 
originate from one kind of oscillations in the transition layer. 

(5) . HBO signals are present in three data sets (30035-01- 
07-00, 93082-01-02-05 and 30035-01-01-00), but they display 
distinct energy dependence. The HBO signal is only detected 
in the original 36-63 channels for 30035-01-07-00 and in the 
original 0-16 channels for 30035-01-01-00. For 93082-01-02- 

05, the HBO is detected in the original 0-63 channels, but not 
in any individual energy band defined in Table [2] It seems that 
the energy dependence of HBO signals changes with the state 
of the source in Scorpius X- 1 . Further investigation on this phe- 
nomenon will be very rewarding but it requires data of suitable 
timing and spetcral resulotion and of good enough statistics. 
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Energy (KeV) 



Fig. 2. The centroid frequency versus the photon energy for 
NBOs in Scorpius X-l. The solid lines are the best fits with 
a two-power-law function to the data. 
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